NASA/TM-2007-2 14623 



Aeroelastic Model Structure Computation for 
Envelope Expansion 


Sunil L. Kukreja 

NASA Dryden Flight Research Center 
Edwards , California 


July 2007 


NASA STI Program ... in Profile 

Since its founding, NASA has been • CONFERENCE PUBLICATION. Collected 

dedicated to the advancement of aeronautics and papers from scientific and technical 


space science. The NASA scientific and technical 
information (STI) program plays a key part in 
helping NASA maintain this important role. 

The NASA STI program is operated under 
the auspices of the Agency Chief Information 
Officer. It collects, organizes, provides for 
archiving, and disseminates NASA’s STI. The 
NASA STI program provides access to the NASA 
Aeronautics and Space Database and its public 
interface, the NASA Technical Report Server, 
thus providing one of the largest collections of 
aeronautical and space science STI in the world. 
Results are published in both non-NASA channels 
and by NASA in the NASA STI Report Series, 
which includes the following report types: 

• TECHNICAL PUBLICATION. Reports 

of completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data 
and information deemed to be of continuing 
reference value. NASA counterpart of peer- 
reviewed formal professional papers but has 
less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary 
or of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to 
NASA’s mission. 

Specialized services also include creating 

custom thesauri, building customized databases, 

and organizing and publishing research results. 

For more information about the NASA 

STI program, see the following: 

• Access the NASA STI program home page 
at http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help @ sti.nasa. gov 

• Fax your question to the NASA STI Help 
Desk at (301) 621-0134 

• Phone the NASA STI Help Desk at 
(301) 621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7115 Standard Drive 

Hanover, MD 21076-1320 


NASA/TM-2007-2 14623 



Aeroelastic Model Structure Computation for 
Envelope Expansion 


Sunil L. Kukreja 

NASA Dryden Flight Research Center 
Edwards , California 


National Aeronautics and 
Space Administration 

Dryden Flight Research Center 
Edwards, California 93523-0273 


July 2007 


NOTICE 

Use of trade names or names of manufacturers in this document does not constitute an official endorsement of such 
products or manufacturers, either expressed or implied, by the National Aeronautics and Space Administration. 


Available from: 

NASA Center for AeroSpace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 
(301) 621-0390 


ABSTRACT 


Structure detection is a procedure for selecting a subset of candidate terms, from a full model 
description, that best describes the observed output. This is a necessary procedure to compute 
an efficient system description which may afford greater insight into the functionality of the 
system or a simpler controller design. Structure computation as a tool for black-box modeling 
may be of critical importance in the development of robust, parsimonious models for the flight-test 
community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion 
that may save significant development time and costs. In this study, a least absolute shrinkage and 
selection operator (LASSO) technique is investigated for computing efficient model descriptions 
of non-linear aeroelastic systems. The LASSO minimises the residual sum of squares with the 
addition of an ^ penalty term on the parameter vector of the traditional f 2 minimisation problem. 
Its use for structure detection is a natural extension of this constrained minimisation approach to 
pseudolinear regression problems which produces some model parameters that are exactly zero 
and, therefore, yields a parsimonious system description. Applicability of this technique for model 
structure computation for the F/A-18 (McDonnell Douglas, now The Boeing Company, Chicago, 
Illinois) Active Aeroelastic Wing project using flight test data is shown for several flight conditions 
(Mach numbers) by identifying a parsimonious system description with a high percent fit for cross- 
validated data. 


NOMENCLATURE 


AAW 

LASSO 

LTI 

NARMAX 


Active Aeroelastic Wing 

least absolute shrinkage and selection operator 

linear time-invariant 

Nonlinear AutoRegressive, Moving Average exogenous 


INTRODUCTION 

System identification, or black-box modeling, is a critical step in aircraft development, 
analysis, and validation for flight worthiness. The development and testing of aircraft typically 
takes many years and requires a considerable expenditure of limited resources. One reason for 
lengthy development time and high costs is inadequate knowledge of an appropriate model type or 
structure to use for parameter estimation. Selection of an insufficient model structure may lead to 
difficulties in parameter estimation, giving estimates with significant biases and/or large variances 
(ref. 1). This often complicates control synthesis or renders it infeasible. The power of using 
structure detection techniques as a tool for model development (i.e. black-box modeling) is that 
it can provide a parsimonious system description that can describe complex aeroelastic behaviour 
over a large operating range. Consequently, this provides models that can be more robust and, 
therefore, reduce development time. 


Moreover, when studying aeroelastic systems it may not be practical to assume that the exact 
model structure is well known a priori. In aerospace systems analysis, one of the main objectives is 
not only to estimate system parameters, but also to gain insight into the structure of the underlying 
system. Therefore, structure computation is of significant relevance and importance to modelling 
and design of aircraft and aerospace vehicles. Structure computation may indicate deficiencies in an 
analytical model and could lead to improved modelling strategies and also provide a parsimonious, 
black-box system description for control synthesis (ref. 2). 

NARMAX MODEL FORM 

Recently, Kukreja and Brenner (ref. 3) showed that NARMAX (Nonlinear AutoRegressive, 
Moving Average exogenous) identification is well-suited to describing aeroelastic phenomena. 
The NARMAX structure is a general parametric form for modeling non-linear systems (ref. 4). 
This structure describes both the stochastic and deterministic components of non-linear systems. 
Many non-linear systems are a special case of the general NARMAX structure (ref. 5). In this 
report, we focus on a special class of NARMAX models: non-linear polynomial models. In eq. (1), 
the polynomial NARMAX structure models the input-output relationship as a non-linear difference 
equation of the form 

z(n) = f l [z(n- 1), - , z(n-n y ),u(n), - ,u(n - n u ),e(n - 1), - ,e(n- n g )\ + e{n) (1) 


/denotes a non-linear mapping, l is the order of the non-linearity, u is the controlled or exogenous 
input, z is the measured output, and e is the uncontrolled input or innovation. This non-linear 
mapping may include a variety of non-linear terms, such as terms raised to an integer power, 
products of present and past inputs, past outputs, or cross-terms. In general, the non-linear mapping, 
/ can be described by a wide variety of non-linear functions such as sigmoids or splines (refs. 5 
and 6). This system description encompasses many forms of non-linear difference equations that 
are linear-in-the-parameters. 

Identifying a NARMAX model requires two things: (1) structure detection and (2) parameter 
estimation. Structure detection can be divided into: (la) model order selection and (lb) selecting 
which parameters to include in the model. We consider model order selection as part of structure 
detection since, theoretically, there are an infinite number of candidate terms that could be 
considered initially. Establishing the model order, then, limits the choice of terms to be considered. 
For the NARMAX model, the system order is defined to be an ordered tuple in eq. (2) as 

0 = K n z n e ( 2 ) 


where n u is the maximum lag on the input, n the maximum lag on the output, n e the maximum lag 
on the error, and l is the maximum non-linearity order. Note that for nonpolynomial NARMAX 
models, l may be simply replaced by a non-linear mapping of some specified class. In this paper, 
we assume that the system order is known. 


2 


STRUCTURE DETECTION 


The structure detection problem is that of selecting the subset of candidate terms that 
best describes the output. Therefore, the parametrisation of a system is still further reduced by 
determining which of the components are required. The maximum number of terms in a NARMAX 
model with n , n u , and n e dynamic terms and /th order non-linearity is defined in eq. (3) as: 


-i P ,+.*,=■ 


i=l 


(3) 


As a result, the number of candidate terms becomes very large for even moderately complex 
models, making structure detection difficult. We define the maximum number of terms, p , as the 
number of candidate terms to be initially considered for identification. Because of the excessive 
parameterisation (the curse of dimensionality), the structure detection problem often leads to 
computationally intractable combinatorial optimisation problems. 

LEAST ABSOLUTE SHRINKAGE AND SELECTION OPERATOR 

Recently, Kukreja et al. (ref. 7) developed a novel approach to this computationally intractable 
combinatorial optimisation problem. The least absolute shrinkage and selection operator (LASSO) 
(ref. 8) is a least-squares-like problem with the addition of an penalty on the parameter vector 
seen in eq. (4) as 


min-||(Z- €>0)||2 + A ||0||j (4) 

e 2 


where II ' II 2 denotes the f 2 -norm, II ' Hi denotes the £j -norm, Z e M ;V x 1 is a vector of outputs, 
0 e K Nxp is a matrix of regressors and 0 e M 2>xl is a vector of unknown system coefficients. 
Note that the regression matrix is a function of the measured input-outputs and unmeasured noise, 
which makes this a pseudolinear regression problem since 0 is (partly) unknown and must be 
estimated along with the parameters. 

The regularisation parameter R. 3 A = controls the trade-off between 

approximation error and sparseness. The LASSO shrinks the least-squares estimator [right-hand 
side eq. (4)] towards 0 and potentially sets Oj = 0 for some j. Consequently, LASSO behaves as a 
structure selection instrument. 

Structure detection can provide useful process insights that can be used in subsequent 
development or refinement of physical models. Therefore, in the sequel, we investigate the 
applicability of LASSO to experimental aircraft data. Specifically, LASSO is used as a structure 
detection tool to assess whether the (i) underlying data is best described by a linear time-invariant 
(LTI) or non-linear model, and (ii) model structure is invariant during envelope expansion. 
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EXPERIMENTAL AIRCRAFT DATA 


The LASSO technique was assessed on experimental flight test data from the F/A-18 Active 
Aeroelastic Wing (AAW) airplane project at the NASA Dryden Flight Research Center. The data 
analysed for this study used collective aileron position input and structural accelerometer response 
output. 


Procedures 

Flight data was gathered during subsonic flutter clearance of the F/A-18 AAW airplane. 
At each flight condition, the aircraft was subjected to band-limited white inputs, with uniform 
distribution and zero-mean. The inputs correspond to collective and differential aileron, collective 
and differential outboard leading edge flap, rudder, and collective stabilator excitations in the 
range of ±0.5 rad and 30 Hz bandwidth for 26 s. 

This report considers accelerometer data measured during the collective aileron sweeps at 
Mach 0.85 and 0.95 both at an altitude of 4,572 m (15,000 ft). The output was taken as the response 
of an accelerometer mounted near the wing leading edge just outside the wing fold (fig. 1). Data 
was sampled at 400 Hz. For analysis, the recorded flight test data was decimated by a factor of 4, 
resulting in a final sampling rate of 100 Hz. 
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Equation (5) identifies the model posed for structure computation was an additive non-linear 
model of the form 


q r 

z(n) = £ 6 v (p(n ) + d w f((p{n )) + + r = p (5) 

V=1 W=1 


The regressors posed for this identification problem were of the form given in eq. (1), up to lag 
order 4, and the non-linear mapping/))), was chosen to be a hyperbolic tangent function (tanh(-)); 
O = [4 4 4 tanh] . It has been observed that aeroelastic structures present a well-defined low frequency 
elastic mode composition (refs. 9 and 10). Therefore, the modal components are represented using 
a fourth-order dynamic expansion (ref. 11). Typically, an assumption is made that a structural non- 
linearity exists in an aeroelastic system that affects not only flutter speed, but also the characteristics 
of the motion (refs. 12 and 13). An example of an aeroelastic non-linearity is a preloaded spring 
tab system. It has been proposed that this type of non-linearity be modelled as a saturator (ref. 13). 
In this paper, we chose to model this type of non-linearity as a hyperbolic tangent, tanh(-), because 
a wing section response typically saturates smoothly. This gave a full model description with 
27 candidate terms. 

The system was identified using scaled hyperbolic tangent functions because the input 
amplitude is less than ± 1 (see fig. 2). The scale factors used for the input, output, and error signals 
were in the range of v = [0.1 1.0] and increased in increments of 0.1. A scaled hyperbolic tangent 
is denoted as tanh(-,v). Models with every possible combination of scale factors were considered 
(i.e. structure computation was performed on 1,000 models). The model which yielded the highest 
cross-validation percent fit was deemed the best-fit model. 

The system was identified applying the LASSO approach. For estimation, A = 2,570 points 
were used from accelerometer response measurements on the left wing. For cross-validation, 
N = 2,570 points were used from data collected at a similar location on the right wing. 
In both the estimation and cross-validation sets, the input was the same collective aileron 
position. The min-max regularisation parameter levels were set to A - n = -10 and A = 0.5 

( 1 (/'mm < X < 1 Q^max ) w i t h a discretisation grid of 1,000 logarithmically spaced A ’s. 

Results 

The results of identifying the AAW data are presented. Figure 2 shows the input-output trials 
used for this analysis. The data represents collective aileron position sequence and structural 
accelerometer response (left wing) used to compute the system structure. 
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(a) Mach 0.85, Alt. 4,572 m (15,000 ft). 
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(b) Mach 0.95, Alt. 4,572 m (15,000 ft). 
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Figure 2. Estimation data. Upper panel: Recorded collective aileron position. Lower panel: 
Recorded structural accelerometer response. 


Equations (6) and (7) depict the model structure computed by the LASSO method 

z(n) = 9 Q + 9^ii(n - 1) + djiiin - 4) + 0 3 tanh(w(?7 - 1),0.4) + § 4 tanh( u{n - 2), 0.4) (6) 

+ tanh(u(n - 3), 0.4) + 9^z(n - 1) + 9jz(n - 2) + 0 g z(?; - 4) + § 9 tanh(z(« - 1),0.4) 

+ 0 1O tanh(z(« - 3), 0.4) + 9 ) ^£(n - 1) + 0 p £(« - 2) + 9^£{n - 4) + 0 J4 tanh(£(n - 1),0.4) 

+ tanh(£(7; - 3), 0.4) 

z(n) = - 1) + / & 2 u(n - 4) + d 3 tanh(w(« - 1),0.4) + d 4 tanh(«(« - 2), 0.4) + 9^z(n - 1) (7) 

+ t7 6 z(/7 - 2) + 9jz(n - 4) + ?7 g tanh(z(« - 1), 0.4) + ?? 9 tanli(z(/7 - 3), 0.4) + d 1Q tanh(z(/7 - 4), 0.4) 

+ d n £(« - 1) + dp£(n - 2) + d 13 £(« - 4) + d |4 tanh(£(« - 1),0.4) + t? 15 tanh(£(« - 3), 0.4) 

+ t 7 16 tanh(£(« - 4), 0.4) 


Equation (6) represents the computed model structure for flight condition Mach 0.85 whereas 
eq. (7) represents the computed model structure for flight condition Mach 0.95. The computed 
model structures are represented by a combination of linear and non-linear, lagged input-output 
terms and contain 15 and 16 terms for Mach 0.85 and 0.95, respectively. Hence, the LASSO 
technique successfully produced a parsimonious model description from the full set of 27 candidate 
terms. 

Figure 3 shows the predicted output for a cross-validation data set for the identified structure 
[(a): eq. (6), (b): eq. (7)]. The upper panel displays the full 26 s time history of the accelerometer 
response recorded on the right wing. The lower panel displays a 10 (Mach 0.85) and 11 (Mach 0.95) 
second slice of the predicted output superimposed on top of the measured output. For Mach 0.85 
[fig. 3(a)] the predicted output accounts for over 95 percent of the measured outputs variance whereas 
for Mach 0.95 [fig. 3(b)], the predicted output accounts for over 98 percent of the measured outputs 
variance. The results demonstrate that the computed model structures are capable of reproducing 
the measured output with high accuracy. 
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(a) Mach 0.85, Alt. 4,572 m (15,000 ft). 
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(b) Mach 0.95, Alt. 4,572 m (15,000 ft). 


8 


Figure 3. Cross-validation data. Upper panels: Full time history of structural accelerometer 
response. Lower panels: Predicted accelerometer response of right wing superimposed on top of 
measured accelerometer output. 


DISCUSSION 


Experimental results demonstrate that structure computation as a tool for black-box modeling 
may be a useful tool for the analysis of dynamic aircraft data. The LASSO successfully reduced the 
number of regressors posed to aircraft aeroelastic data, yielding a parsimonious model structure 
for each data set. Additionally, these parsimonious structures were capable of predicting a large 
portion of the cross-validation data collected on the adjacent wing and with a different sensor. This 
suggests that the identified structures and parameters explain the data well. Using percent fit alone 
as an indicator of model goodness could lead to incorrect interpretations of model validity. In many 
cases, however, for non-linear models this may be the only indicator that is readily available. 

In this work, the results show that whereas the linear dynamics remained invariant for 
both flight conditions available for analysis, the non-linear dynamics changed as Mach number 
increased. For Mach 0.85 the model [eq. (6)] displayed a richer non-linear dynamic description 
associated with the input, possibly a result of actuation stiction. When the Mach number was 
increased to 0.95, with other flight conditions held invariant, the model [eq. (7)] demonstrated 
a richer non-linear dynamic description associated with the output, possibly associated with the 
ailerons. It is difficult to make definitive comments on the underlying physics responsible for this 
behaviour without extensive analysis of different flight conditions. The important points to note 
are, this study suggests (i) non-linear models are appropriate to describe the dynamics behaviour 
of advanced aircraft, and (ii) models describing aircraft dynamics vary with flight condition. This 
suggests non-linear modelling may afford a robust and parsimonious system description over a 
larger operating regime and models used for prediction (e.g. control) are not invariant for all 
flight conditions. This may hold significant implications for aircraft development, saving time 
and money. For this study, only hyperbolic tangent with fourth-order input, output and error lag 
was used as a basis function to explain the non-linear behaviour of the F/A-18 AAW data. Clearly, 
different basis functions and a higher dynamic order (lag-order) should be investigated to determine 
if another basis could produce accurate model predictions with reduced complexity. Moreover, 
further studies are necessary to evaluate whether the model structure is invariant under different 
operating conditions, such as altitude and model parameterisations. 

This study illustrates the usefulness of structure detection as an approach to compute a 
parsimonious model of a highly complex non-linear process, as demonstrated with experimental 
data of aircraft aeroelastic dynamics. Moreover, analysis of input-output data can provide useful 
process insights that can be used in subsequent development or refinement of physical models. In 
particular, morphological models are based on assumptions (e.g. these effects are important and 
those are negligible) which may be incorrect (refs. 14 and 15). A structure computation approach 
to model identification may help uncover such surprises. 
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CONCLUSIONS 


The LASSO is a novel approach for detecting the structure of over-parameterised non-linear 
models. These results may have practical significance in the analysis of aircraft dynamics during 
envelope expansion and could lead to more efficient control strategies. In addition, this technique 
could allow greater insight into the functionality of various systems dynamics, by providing an 
easily interpretable quantitative model. 
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